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j_( . Abstract. By adding a weak magnetic field to a spherically symmetric fluid configuration that 

' caricatures a stalled shock in the post-bounce supernova environment, we explore the capacity of 

C3 , the stationary accretion shock instability (SASI) to generate magnetic fields. The SASI develops 

upon perturbation of the initial condition, and the ensuing flow generates — in the absence of 
rotation — dynamically significant magnetic fields (~ 10 15 G) on a time scale that is relevant 
for the explosion mechanism of core-collapse supernovae. We describe our model, present some 
recent results, and discuss their potential relevance for supernova models. 

o ' 



O ■ 1. Introduction 

Magnetic fields in the context of core-collapse supernovae — the final evolutionary stage of massive 
stars — have recently attracted considerable interest (see [13 [7], and references therein) due, in 
part, to the generically observed asphericity of supernova explosions [16] , 

In particular, simulations of rapidly rotating progenitors followed for several hundred 
!— i ' milliseconds past core bounce demonstrate the ability of magnetic fields to drive and collimate 

bipolar explosions [7j: Differential rotation near the surface of the rapidly rotating proto-neutron 
star (PNS) winds and organizes poloidal magnetic fields into 'magnetic towers' along the rotation 
axis, consisting mainly of toroidal fields. Although this mechanism is compelling, current models 
rely on rotation rates — and an initial field strength and topology — that are incompatible with 
stellar evolution calculations that take magnetic effects into account [10j . 

It is often assumed that the magneto-rotational instability (MRI), which may well operate 
in core-collapse supernovae [2], will bridge the gap and make the winding mechanism work. 
Capturing the MRI in numerical simulations requires higher spatial resolution than has been 
possible in realistic supernova simulations. The MRI has been observed in a simplified model of 
stellar collapse [15], but the potential impact of the MRI (and magnetic fields generally) on the 
core-collapse supernova explosion mechanism remains unclear. 

A different mechanism that may account for aspherical supernovae is the SASI, a 
hydro dynamical instability of the stalled supernova shock wave. Idealized two-dimensional (2D) 
and three-dimensional (3D) simulations have shown the SASI to have a dramatic impact on the 
postshock flow [31[4]: The £ = 1 (sloshing) mode is seen to dominate 2D simulations, while the 
m = 1 (spiral) mode dominates in 3D and may generate significant pulsar spin. The SASI is 
also seen to play an important role in the dynamics of more realistic models [3 Q3] • 



Given its impact on the postshock flow, it is natural to explore the influence of the SASI on 
magnetic fields (and vice-versa) in the supernova environment. By extending an axisymmetric 
model j3] to include magnetic fields we have initiated an investigation of the evolution of 
magnetic fields when driven by the SASI: Starting from a spherically symmetric, non-rotating 
configuration resembling a stalled supernova shock wave, we find that axisymmetric SASI- 
induced flows are able to amplify magnetic fields to dynamically significant levels within the 
time frame a supernova is expected to reach explosive conditions. 

2. Physical model and numerical solution 

In order to study magnetic field generation by the SASI, and the potential relevance to the 
core-collapse supernova explosion, we adopt a highly idealized description of the post-bounce 
supernova environment [3]: A steady-state, spherically symmetric accretion shock is located at 
r = R Sh = 200 km from the center of the PNS. Matter at larger radii is falling into the shock at 
the free-fall speed, u = \/2GM/r, where G is Newton's constant, and the mass M of the central 
object is set to 1.2 M@. A fixed accretion rate of 0.36 Mq s _1 is employed, and a constant, 
highly supersonic Mach number of 300 is used to set the pressure in the pre-shock gas. The 
gravitational potential is given by the point-mass formula <3? = —GM/r. The Rankine-Hugoniot 
conditions determine the hydrodynamic state just inside the shock, and the Bernoulli equation 
is solved for the structure from the shock to the surface of the PNS, which serves as the inner 
boundary and is located at r = R PNS = 40 km. The initial magnetic field is that of a split 
monopole: B = B r e r , where B r = sign (cos 9) x Bq (r/R PNS ) 2 , and 9 is the polar angle in the 
spherical coordinate system. The strength of the magnetic field at r = R PNS , Bq, is set to 
10 10 G in our base model. The magnetic field lines are initially parallel to the accreting flow and 
the configuration is consistent with the steady state initial condition. The system is described 
by the non-relativistic ideal magnetohydro dynamic (MHD) equations, which are closed with 
a polytropic equation of state (EoS) with the ratio of specific heats 7 set to 4/3, similar to 
collapsing matter in a supernova. 

We have implemented a time-explicit, second-order, semi-discrete central- upwind scheme [11] 
in our code, GenASiS, for the integration of the time-dependent ideal MHD equations. The 
equations are formulated in their integral form, appropriate for finite-volume shock-capturing 
methods. We use the HLL formulae [13] to compute the numerical fluxes and electric fields, 
located on faces and edges of a computational cell, respectively. Second-order spatial accuracy 
in smooth regions, while maintaining a non-oscillatory behavior near shocks and discontinuities, 
is achieved by linear interpolation subject to a generalized minmod limiter [12]. The method 
of constrained transport [9] is used to maintain the divergence-free condition on the magnetic 
field during the evolution of the MHD equations. Second-order temporal accuracy is obtained 
with a two-step Runge-Kutta time step algorithm [12]. The equations are solved on a uniform 
Cartesian grid (Ax = Ay = Az = Al), and in order to solve axially symmetric problems we 
have extended the 'Cartoon' method [1] to the MHD case. 

The results presented in this study have been computed with a spatial resolution given by 
Al/R 

pns — 0.06. At this resolution we obtain solutions that are consistent with the MHD 
equations in the sense that the total energy in the system is conserved at an acceptable level, 
and that the unperturbed model can indefinitely maintain (modulo a slight initial relaxation) 
the steady state of the initial condition, without triggering the SASI from numerical noise. 

3. Results 

In order to initiate the SASI we perturb the initial condition by placing two dense rings in 
the material ahead of the shock. As the density enhancements pass through the shock they 
introduce non-radial perturbations that grow and initiate the SASI. The £ = 1 mode quickly 
manifests itself as an up-down 'sloshing' of the global accretion shock configuration. In Figure 



(HJ) we plot a time series illustrating the evolution of the SASI from the onset of the calculation 
and well into the non-linear phase: We show a color plot of the magnetic field strength, |B|, at 
four different times during the evolution. The snapshots shown are taken at 294 ms (a), 336 ms 
(b), 424 ms (c), and 650 ms (d). 



Figure 1. Color plot of the 
magnitude of the magnetic field, 
|B|, during the evolution of the 
SASI. Panel (a) is sampled at 
294 ms, panel (b) at 336 ms, 
while lower panels, (c) and (d), 
are sampled at 424 ms and 
650 ms, respectively. The color- 
bar in the upper right (valid 
for all panels) indicates the 
magnetic field strength on a log 
scale, extending from 10 6 G to 
10 15 G. Thin lines are contours 
of constant density; starting 
with the innermost contour we 
plot density contours at p = 
10 10 , 10 9 , 3 x 10 8 , and 6 x 10 7 
g cm -3 . The thick contour 
indicates the location of the 
shock surface. Notice the factor 
of two difference in the spatial 
scale used in the upper and 
lower panels of the plot. 

During SASI development the shock starts to deviate from spherical symmetry, the accreting 
flow hits the shock at an oblique angle, and significant non-radial flows are introduced inside 
the shocked cavity. The lateral flow carries magnetic field toward the symmetry axis, where the 
strength of the field increases due to compression: At 294 ms the strongest field is found around 
the north pole, with -B max ~ 2 x 10 12 G. At the south pole the field strength is about 3 x 10 11 G. 
Some 40 ms later, at 336 ms, the situation is reversed, and the strongest field, -B max ~ 2 x 10 12 G, 
is found around the south pole. 

An internal shock forms in the postshock gas around 360 ms, and is connected to the accretion 
shock at a triple point [1] (visible as a kink in the shock surface). The internal shock propagates 
in the north-south direction in panel (c). Ahead of the internal shock, accretion-powered flow 
plunges toward the PNS and the axis of symmetry, causing further field growth. At 424 ms 
columns of stronger magnetic field are seen in both polar regions, with maxima of some 10 G. 
At 650 ms extended regions of strong magnetic field have formed around both the north and 
south poles. The strength of the magnetic field in these regions is above 10 15 G near the surface 
of the PNS, and stays above 10 14 G out to about 100 km beyond the surface of the PNS. The 
magnetic field is these regions is strong enough to influence the dynamics: From the density 
contours in panel (c) we see that low-density funnels have formed along the polar axis. 

The impact of the SASI- induced flow on the magnetic field is further illustrated in Figure (J2]) , 
in which we plot subsets of the integrated kinetic (red lines) and magnetic (black lines) energies 
inside the shocked cavity. The solid red line represents the total kinetic energy inside the shock, 





Figure 2. Kinetic (red curves) and magnetic 
(black curves) energies (see text for detail): 

Elt ( ), EL ( ), (-- •-), and 

EL ( )• *£«( ),^l g (---"), and 

EL* ( )• We also plot for an 

unperturbed model ( — • — ). (1 B=10 51 erg.) 



Figure 3. Magnetic energy inside the 
shocked cavity, E^ , for models with various 
values of initial magnetic field strength Bq 

at the PNS surface: 10 8 G ( ), 10 9 G 

( ), 10 10 G ( ), and 10 11 G ( ), 

respectively. 



i^kin = J v u • udV, where V Sh is the total volume encompassed by the shock. The red dash- 
dot line represents the kinetic energy of the lateral motion, E^ in = f v \pu 2 e dV. The dashed 
and dotted red lines represent the kinetic energy around the northern and southern hemispheres, 
respectively: E^ in = J v ^ ^pu-udV and E^ in = J Vs ^pu-udV. Here V N(S) is the volume bounded 
by the surface of the PNS, the surface of the accretion shock, and the surface a cylinder with 
radius -R PN s, centered on the polar axis in the northern (southern) hemisphere. The solid black 
line represents the magnetic energy inside the accretion shock, E^ g = j Vi ■ BdV, while 

the dashed and dotted lines represent the magnetic energy in cylinders around the north and 
south pole: E^ a? = J Vn • B dV and E^ = J Vs ^B • B dV respectively where p is the 

vacuum magnetic permeability. For comparison, we also plot the magnetic energy inside the 
shock, E^ , for an unperturbed model (black dash-dot line). 

The non-radial perturbations lead to the growth of the kinetic energy of the lateral motion 
inside the shock. This is seen in the rise in E® in for t < 400 ms. The sloshing nature of the 
I = 1 mode is also revealed in the alternating pattern of the kinetic energies contained in the 
cylinders centered around the northern and southern hemispheres. Note the correspondence 
of rise in magnetic energies with decline in kinetic energies. This persists to about 450 ms, 
when the flow takes on the character of the I = 2 mode. From the evolution of the magnetic 
energy we can see that the SASI has a quite dramatic impact on the magnetic field. The early 
evolution is characterized by relatively modest oscillations in the strength of the magnetic field 
superimposed on a trend of gradual overall increase, but starting at about 360 ms, around the 
formation of the internal shock, the magnetic energy increases by six orders of magnitude in 
roughly 100 ms. Beyond 450 ms the magnetic energy saturates at a level somewhat below the 
total kinetic energy inside the shock. The magnetic field is strongly aligned with the symmetry 
axis and the magnetic energy is in approximate local equipartition with the kinetic energy in 



the polar regions. 

We have also done calculations with different initial field strengths. In Figure ([3|) we plot the 
magnetic energy inside the accretion shock, E^ 1 , for models where we have varied Bq: 10 8 G 
(dashed line), 10 9 G (dotted line), 10 10 G (solid line, the model shown in Figures ([T]) and fl2J)), 
and 10 11 G (dash-dot line). These calculations show that our result is insensitive to the strength 
of the initial field: The magnetic energies saturate at the same level for all models. 



4. Summary and discussion 

By considering an idealized model of a stalled supernova shock wave we have found that the 
SASI-induced flow beneath the shock is able to generate dynamically significant magnetic fields, 
in excess of 10 15 G, on a time scale relevant to the supernova explosion, independent of the 
initial field strength. The influence of the SASI on the evolution of magnetic fields in core- 
collapse supernovae will need to be investigated in more realistic models, but our idealized 
models offer a proof of principle that (1) magnetic fields may be amplified to significant strengths 
even in the absence of rotation and (2) compression may serve, at least in axisymmetry, as an 
effective mechanism for rapid field amplification to dynamically significant levels, contrary to 
conclusions reached in past MHD collapse calculations. Moreover, current 2D models [5| l6j 18} [T4] 
demonstrate that the SASI is common to core-collapse supernova models. Our results then 
suggest that SASI magnetic field amplification and, in particular, the formation of dynamically 
significant fields may be common as well, and not confined to special cases — for example, to cases 
of rapid progenitor rotation. The key question now is: What happens in three dimensions? The 
amplification we document here relies on the continual focusing of fluid flow toward the polar 
axis, in part an artifact of our assumed axisymmetry. Hydrodynamics simulations of the SASI 
in 3D yield a significantly different and more complex postshock flow [4], whose consequences 
for stellar core magnetic field amplification remain to be investigated. 
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